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Abstract. - We study the distribution of overlaps of glassy minima, taking proper care of 
residual symmetries of the system. Ensembles of locally stable, low lying glassy states are 
efficiently generated by rapid cooling from the liquid phase which has been equilibrated at a 
temperature T Iun . Varying T Iun , we observe a transition from a regime where a broad range 
of states are sampled to a regime where the system is almost always trapped in a metastable 
glassy state. We do not observe any structure in the distribution of overlaps of glassy minima, 
but find only very weak correlations, comparable in size to those of two liquid configurations. 



The phenomenology of supercooled liquids, the glass transition and the glassy phase have 
been related to results of analytical calculations for microscopic models of spin glasses JIJ || . 
The latter bear no obvious relevance to structural glasses, hence, it is particularly important to 
see, whether the analogy can be supported quantitatively. One of the most prominent features 
of those spin-glass models, which have been put forward as models of structural glasses, is 
one step replica symmetry breaking, resulting in a bimodal distribution of overlaps between 
pure states. Generally speaking, one would expect to see non-trivial correlations among glassy 
states, if mean field theory of spin glasses were to apply. 

The liquid to glass transition as it is observed e.g. in metallic glasses, is a non-equilibrium 
phenomenon with the crystalline state being the true ground state. Glassy states are in general 
metastable states, sometimes called quasi-ergodic components, which are mutually inaccessible 
on experimental time scales, while equilibrium within a quasi-ergodic component is reached 
on much shorter time scales. In contrast to mean field models of spin glasses the ensemble 
of metastable states depends on the preparation method and the appropriate weights for 
metastable states are not unique. The assumption of an idealised preparation method, which 
cools the liquid arbitrarily slowly and at the same time prevents crystallization, has never 
been verified and seems contradictory to us. We therefore prefer to use a different idealised 
preparation ensemble which corresponds to infinitely rapid cooling from liquid equilibrium 
states. This method was first introduced by Stillinger and Weber |3, 
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Distance Measure in Configuration Space. - We study correlations among an ensemble 
of such metastable states, labelled by a. A glassy state a is characterized by a spontaneous 
breaking of translational invariance as signalled by a non-zero expectation value of the Fourier- 
components (exp (ikfi)) a 7^ of the local density. Here fi denotes the position of the i-th 
particle, i = 1, . . . ,N. The local density is analogous to the local magnetization, which is 
non-zero in the spin glass phase, (s,) a ^ 0, indicating the spontaneous breaking of the Ising 
symmetry. An overlap between two glassy configurations can be defined as 

1 N 

Q a -\k) = - ]T(c X p (ikn)) a (cx P {-ikn)) b (1) 
1=1 

in naive analogy to the spin glass order parameter Q ab = jjJ2iLi( s i)a( s i)b- For the glassy 
state, however, Q a ' b (k) is not a good measure, because it depends on the labelling of the 
particles. In particular two identical configurations with different labels can have a very small 
overlap according to the above definition. Permutations of identical particles do not give rise to 
new states in the space of physically distinguishable configurations. Hence, we need a distance 
measure which is invariant under all permutations ir of identical particles in one state only. 
This can be achieved by maximizing over all such permutations 

1 N 

Q a - b (k) = max — ^(cxp {ikr t )) a (exp (-ikf v{i) )) h . (2) 

i=l 

This overlap has the required invariance property, but its computation requires the enumer- 
ation of all permutations, which is computationally hard and can at most be done for very 
small systems. 

The computationally hard problem can be avoided by working with densities. We specialize 
to zero temperature and only consider configurations of locally minimal energy. Such a 
configuration is uniquely characterized by a point measure p a (r) — (^/N)^2f =1 S(r — rf), 
which by construction is independent of the labelling of the particles. To obtain a distance 
measure in configuration space, we regularise the point particle densities 5(f) to homogeneous 
77-spheres 

5(f) ->■ A r) (r) = - |f|), (3) 

Here V v — Airrj 3 /3 is the volume of the 77-sphere and we choose 1r\ < r^ in = mini.j(|r^ a — rj^l) 
to guarantee that the spheres do not overlap and that the positions ff- can still be uniquely 
reconstructed from the regularised density. Such a regularisation is necessary because products 
of 5 functions are ill defined and the probability of coincidence of any two 3-dimensional vectors 
which are distributed smoothly in space is equal to zero. 

A natural distance measure between two regularised densities is 

q^ = ^ r d z x ST A v (x-xnA v (x-x b ). (4) 

If we restrict r] to even smaller values, 47/ < r min :— min(r^ in , r^ in ), then each sphere in 
state a can overlap with at most one sphere in state b. Consequently the product of two local 
regularised densities is only non-zero, for the permutation j — 7r(i), which identifies particles 
within the same ^-sphere. The overlap is then given by the volume fraction of overlapping 
spheres in the two configurations, so that < Q a ' b < 1. 



Energy Landscape of LJ Glasses 



3 



It is obvious that 77 cannot be made arbitrarily small, because then the distance measure 
looses its ability to discriminate between different structures. Very small 77 will simply lead 
to maximal distances in most cases and ultimately reproduce the point measure. We have 
varied 77 over the range < 77 < rmin/4 and computed the number of overlapping spheres for 
two configurations which are random and two configurations which are similar (about 80% 
overlap) . We find that the number of overlapping spheres is hardly sensitive to the choice of 77 
in the range r m i n /8 < 77 < r nun/4. For 77 < r min /8 the number of overlapping spheres decreases 
drastically with 77 even for the strongly correlated configurations so that Q a b can no longer 
discriminate between strongly and poorly correlated states. In the following we fix 77 = r m - ln /5. 

Other distance measures have been used in the literature [El , in particular 



which requires the solution of a hard computational problem. It has the further disadvantage 
to be not directly related to quantities which may be obtained from experiment (like the 
density) or which appear naturally in any of the existing theories of supercooled liquids and 
structural glasses. 

A system in free space with pairwise interactions possesses translational and rotational 
symmetries, which may be broken explicitly by boundary conditions. The most convenient 
boundary conditions for computer simulations, namely a cubic box with periodic boundary 
conditions, leave a subgroup of the free space symmetries unbroken. This subgroup S is 
generated from the translations and the cubic point group symmetries. As a consequence, there 
exists a whole orbit of symmetry related, degenerate minima of H to every single minimum 
found in the simulations. 

We are interested in the rate of occurrence, i.e. the probability distribution, of overlaps 
P(Q) = J2 a b-PaPb<>{Q — Q a ' b ). Here the summation includes all symmetry related states, 
which all have the same rate of occurrence P a . The large number of symmetry related states 
imposes two severe difficulties. First, many translated states (out of a continuum of states) 
have to be generated to obtain a reliable estimate of P(Q) from numerical simulations. Second, 
the Q values generated by all symmetry related states may scatter considerably in the interval 
< Q < 1. This may smear out structures in P(Q) which would be clearly marked in an 
ensemble without residual symmetries. To avoid this problem we have broken the translational 
symmetry explicitly by fixing the centre-of-mass. Since total momentum is conserved in our 
simulations, this can be achieved by an appropriate choice of initial conditions (see below). The 
cubic point group symmetry has been left untouched: for each minimum which was found by 
our algorithm, we have generated 48 equivalent states and included them in the histogram for 
Qa,b rp Q c^gcj^ whether this procedure smears out relevant structure, we have also maximized 
the overlap over all 48 equivalent states. 

Simulations. - The system under consideration is a binary mixture of large (L) and small 
(S) particles with 80% large and 20% small particles. Small and large particles only differ in 
diameter, but have the same mass. They interact via a Lennard- Jones potential of the form 
U a p{r) = 4,e a p[(a a p/r) 12 — [o a p/r) G ]. All results are given in reduced units, where oll was 
used as the length unit and as the energy unit. The other values of e and a were chosen as 
follows: €ls — 1.5, ctls = 0.8, egg — 0.5, agg — 0.88. The systems were kept at a fixed density 
p w 1.2. Periodic boundary conditions have been applied and the potential has been truncated 
appropriately according to the minimum image rule M , and shifted to zero at the respective 
cutoff. For the systems with N = 30 particles the cut-off is r c = 1.43 and for N = 60 the 
cut-off is at r c = 1.7. The choice of the Lennard- Jones parameters follows recent simulations 
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of Lennard- Jones glasses |9|, |l0[; it is known to suppress recrystallization of the system on 
molecular dynamics time scales. The glass transition is supposed to occur at the temperature 
T g « 0.45 ||. Throughout this study we present results for systems with N = 30 and N = 60 
particles, noticing that most of the results have been verified for N = 100. 

Initially N atoms are placed randomly inside a cubic simulation box of side- length L. The 
following steps are performed repeatedly: 

1 Heat up the system to a temperature T run . 

2 Let the system evolve for a time r run using molecular dynamics. 

3 Locate the nearest local minimum by quenching down the system to T = OK using the 
steepest descent path. 

Explicit breaking of translational invariance can be achieved by fixing the centre-of-mass of 
the system. A unique definition of the centre-of-mass R s is not obvious for periodic boundary 
conditions. We start from the extended zone scheme representation of the simulation box. 
The infinite periodic density pattern is cut off at finite volume V, so that the cubic point 
group symmetry stays intact. The centre-of-mass is the same for all such volumes, including 
the smallest one, the simulation box itself. The regularised R s remains unchanged in the limit 
V — ► oo. 

Each run is initialized with the same positions and different velocities of the particles. 
The latter are drawn from a Maxwell-Boltzmann distribution with temperature T run , subject 
to the constraint that the total momentum vanishes. The total momentum is conserved in 
the molecular dynamics simulation as well as in the steepest descent algorithm. Thereby we 
generate an ensemble, such that in all states the centre-of-mass stays fixed. In the following 
we shall also consider states which are related by cubic point group symmetries. To guarantee 
that the centre of mass remains unaffected by point group operations we choose R s = 0. We 
use standard molecular dynamics with the velocity form of the Verlet algorithm . Steepest 
descent is achieved with the conjugate gradient algorithm ||. 

Results. - Simulating small systems requires careful checks to avoid sampling crystalline 
configurations. First, it can be quite helpful to look at the structures. Second, we have applied 
a common neighbor analysis JllJ to detect the amount of polytetrahedral order (amorphous) 
and closed-packed structures (crystalline). Occasionally the structural stability of low energy 
states has been verified for temperatures below the glass transition. 

In fig. [I we show a histogram of the energy per particle of locally minimal states for N = 30 
and N = 60. For each system size 1000 minima have been generated, starting from T run = 0.5. 
For the small system, crystalline states have the dominating basins of attraction, whereas for 
the larger system we found very few crystalline states (3 out of 1000). 

Energy Landscape. The set of sampled local minima depends on T run . We have generated 
200 configurations for N = 60 particles and fixed T run chosen in the interval 0.3 < T run < 
T = 1.8. For each T run we have computed the mean energy, the lowest energy found and the 
variance of the energy. The results are shown in fig. |2| For each configuration the length of 
the molecular dynamics simulation has been chosen as r run = 0.5 * 10 6 . This is sufficient to 
allow for equilibration at T run > T g . 

The mean energies per particle of the sampled minima are approximately constant ((e) rs 
eh = —5.2) at high temperatures and decrease significantly within a transition region. The 
lowest energy minima e; « —5.45 are found for T run rs T g , whereas for T run < T g sampling 
is restricted to one quasi-ergodic component. Thus extensive sampling of locally stable states 
is most effective for T run > T g . At high temperatures, most of the weight is found in a 



Energy Landscape of LJ Glasses 



5 




Fig. 1. - The energy spectrum of the minima sampled from T run = 0.5 for N = 30 (left) and N = 60 
particles (right) is shown as a histogram for 1000 states for each system size. 

Fig. 2. - Mean energy and variance, as calculated by averaging over 200 configurations at each 
temperature T run . The lower graph in the upper panel is the lowest energy found in the set of 
configurations at each temperature. 



broad band of locally stable states with high energies around en ■ The low energy structures 
are not sampled, although these minima are favoured by a relative Boltzmann factor ~ 
exp (N(eh — ei)/fcsT run ). This implies that the configuration space of the high energy basins 
is much larger than that of the low energy basins, so that it overcompensates the relative 
Boltzmann factor. When T run is lowered to temperatures around T g , the Boltzmann factor 
becomes more important and eventually overruns the configuration space volumes of the eh 
states. These results are in agreement with recent work of Sastry et al. who showed that 
the transition is accompanied by the onset of typical, glass like relaxation phenomena. 

To assess the performance of our algorithm in sampling low lying energy states quanti- 
tatively, we compare the lowest energies sampled with the results of a recent study on the 
optimization properties of traditional simulation algorithms applied to the Lennard-Joncs 
glasses p3fl . The authors reported on a lowest energy configuration found for N — 60 at 
e = —5.38 by cooling down slowly from the liquid phase using standard molecular dynamics. 
The lowest state shown in fig. || is at e = —5.45. The difference in energy 5e — 0.07 corresponds 
to a factor of 100 in computer time according to the estimate given in refs. jl3|, |l(J. For 
Trun < 0.5 around 15% of all states sampled using our method are lower in energy then 
e = -5.38. 

Distribution of overlaps. To investigate the distribution of overlaps we have generated 
2000 amorphous states for the N = 60 particles system keeping the centre-of-mass fixed. 
Each molecular dynamics trajectory has been initialized at T run = 0.5 and has been run for 
r run = 10 6 time steps. We compute the overlap of the large particles for all pairs (a,b). This 
implies that the weights P a are approximated by the rate of occurrence of state a in our 
simulation. Several arguments support our choice of ensemble: 1) The generated states are 
as low in energy as the best ones from other optimization procedures. 2) In agreement with 
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recent simulations of Sastry et al. [pL 2|| , we find a qualitative change in the properties of sampled 
states at a well defined temperature T g « 0.5, such that above T g a broad range of states are 
sampled, whereas below T g the system is landscape dominated, i.e. trapped in a metastable 
state. 3) The Stillinger- Weber method is controlled and reproducible and sufficiently simple 
to allow for analytical calculations. 

A histogram of overlaps is shown in fig. [| We observe a most probable value around 
Qa,b 2 anc | n0 significant structure in the distribution. It is instructive to also compute 
the distribution of overlaps of liquid configurations. This distribution is compared to the 
histogram for glassy states in fig. ^ The distribution of overlaps for glassy states is 35% 
broader and centered at a 20% higher value. Otherwise no significant changes occur. Thus we 
conclude that the glassy minima are as different as they can be. This conclusion also holds, if 
we restrict our ensemble of glassy states to the lowest energy states. 




Fig. 3. - The overlap distribution for all glassy states in comparison to the overlap distribution of 
random states, generated as configurations of the liquid phase at high temperatures. 

Fig. 4. - The distribution of all overlap computed by either maximizing over the 48 symmetry related 
states (solid line) or by including all symmetry related states with the same weight (dotted line). 



The residual cubic point group symmetry does not affect the overlap distribution signifi- 
cantly. We have calculated overlaps by either maximizing over the 48 symmetry related states 
or by including all symmetry related states with the same weight. The resulting distributions 
are shown in Fig. ^. As one would expect the distribution obtained by maximizing over all 
rotations is peaked at a higher value of Q, but does not reveal any additional structure. 

Recently Coluzzi and Parisi || computed the distribution of distances P(D), using the 
distance measure of eq. (^) and simulating small systems (28 to 36 particles) confined by soft 
walls. They find highly non-trivial distributions P{D) in the glassy phase, which strongly 
depend on particle number. We have also applied our method to generate an ensemble of low 
energy states by rapid cooling for the system of ref. 13. A straightforward evaluation of P(Q) 
yields the bimodal distribution, shown in fig. ^| A closer inspection, however, reveals that 
the structure is due to many imperfect crystalline states. Soft walls increase the probability 
for crystalline states, as compared to periodic boundary conditions, and for small systems 




Fig. 5. - The overlap distribution for a small system of N = 30 particles, confined by soft walls. The 
dominant peak at high Q is due to overlaps between two crystalline states, whereas the lower values of 
Q correspond to overlaps between two amorphous states or one amorphous and one crystalline state. 
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